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Abstract 

The volume fluctuations in the steady state reached by a vibrated granular gas of hard particles 
confined by a movable piston on the top are investigated by means of event driven simulations. Also, 
a compressibility factor, measuring the response in volume of the system to a change in the mass 
of the piston, is introduced and measured. From the second moment of the volume fluctuations 
and the compressibility factor, an effective temperature is defined, by using the same relation as 
obeyed by equilibrium molecular systems. The interpretation of this effective temperature and its 
relationship with the granular temperature of the gas, and also with the velocity fluctuations of 
the movable piston, is discussed. It is found that the ratio of the temperature based on the volume 
fluctuations to the temperature based on the piston kinetic energy, obeys simple dependencies on 
the inelasticity and on the piston-particle mass ratio. 

PACS numbers: 45.70.-n,47.70.Nd,51.10.+y 
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I. INTRODUCTION 

Granular matter in general jl, 2], and granular gases in particular js-5], have been recently 
the object of intensive theoretical and experimental research. Not in the least this is be- 
cause of the opportunity they offer to investigate fundamental questions of non-equilibrium 
physics. One of these issues is the nature and properties of non-equilibrium fluctuations and 
their relevance for the description of the macroscopic behavior of the system. The simplest 
fluctuations that can be considered are those of global properties of a system. The total 
energy fluctuations of an isolated granular gas, modeled as an ensemble of smooth inelastic 
hard spheres, have already been investigated elsewhere . They are specific of granular 
systems, and a consequence of the localized character of the energy dissipation in collisions. 

In this paper, the volume fluctuations of a vibrated granular gas are studied by means of 
event driven simulations js, sj. The volume of the system changes because the wall on top 
of the grains is a piston that can move in the vertical direction. The equilibrium position 
of the piston is determined by equating its weight per unit of area with the pressure of the 
granular gas just below it. On the other hand, the properties of the fluctuations of the 
piston around that position, and therefore the volume fluctuations of the system, are not 
known. If the system were a molecular gas at equilibrium, therefore without being energized 
by vibration, the total volume fluctuations are known to be Gaussian, and with a second 
moment that is proportional to the temperature and the isothermal compressibility of the 



gas 



10]. There is no reason to expect the above properties to hold in granular gases, that 



are inherently in non-equilibrium states. Even more, to give a meaning to them, a first 
point to be addressed is to specify the notion of temperature to be used. Actually, the 
equilibrium expression for the volume fluctuations can be used to define a temperature-like 
quantity, that can be experimentally measured. This is somehow an extension of the usual 
way of defining "effective temperatures" from the extension to non-equilibrium states of 
the Fluctuation-Dissipation Theorem. A revision of these ideas in the context of granular 
media is given in Uj. See also 12] for a discussion of the validity of the Einstein relation 
in externally driven granular gases. 

The emerging natural question is whether the effective temperature defined from the 
relationship between volume fluctuations and compressibility has some other conceptual 
interpretation, and if it is related to other possible sensible definitions of temperature in the 
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system. The first obvious candidate to be considered is the granular temperature of the gas, 
defined from the second moment of the velocity distribution of the grains. Actually, this 
parameter is known to play in the hydrodynamic description of granular gases a role similar 
to the usual temperature in molecular hydrodynamics [13j . But there is another temperature 
parameter that is relevant for the description of the system, the one defined from the second 
moment of the velocity distribution of the movable piston. Both temperatures, defined from 
the velocity of the gas and of the piston, respectively, are the same in equilibrium systems, 
but they can differ strongly in granular systems, as a manifestation of the violation of energy 
equipartition ]j| . Clarifying the relationship between the effective temperature defined from 
the volume fluctuations, the granular temperature, and the temperature parameter of the 
piston, is one of the aims of this work. 

The structure of the remaining of this paper is as follows. In Sec. [Til the system is 
described and the macroscopic steady state to be considered is characterized. Also, the 
parameter region to be investigated is specified. It is shown than when the mass of the 
piston is much larger than the mass of the grains, the volume fluctuations of the system 
exhibit a Gaussian distribution. The second moment of this distribution depends much 
stronger on the inelasticity of collisions between particles than on the elastic or inelastic 
character of the collisions of the particles with the piston. In Sec. IHH a compressibility 
factor is defined for the granular gas. As usual, it is a measurement of the change in the 
volume of the system as a consequence of a change in the external pressure under well 
defined constrains. From the values of the second moment of the volume fluctuations and 
of the compressibility factor, an effective temperature is defined, as indicated above. This 
temperature turns out to be quite simply related with the temperature parameter of the 
piston, but its relationship with the granular temperature of the gas appears to be very 
intricate. Finally, Sec. IIVI contains a short summary of the main results in the paper and 
some general comments. 



II. STEADY FLUCTUATIONS OF THE POSITION OF THE PISTON 

Consider a system composed by N smooth inelastic hard disks of mass m and diameter 
d, in presence of gravity, and confined by a movable piston of mass M located on the top. 
By definition, the piston can only move in the direction of the gravity field. Moreover, 
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there is no friction between the piston and the lateral walls of the vessel containing the gas. 
The system of particles is kept ffuidized and at low density by injecting energy through the 
bottom wall, which is vibrating. Inelasticity in collisions between particles is modeled by 
means of a constant coefficient of normal restitution, a, defined in the interval < a < 1. 
Thus when two particles i and j collide their velocities change instantaneously from the 
initial values Vi, vj to the post-collisional ones given by 

1 + a 



v- = Vi 



2 

1 + a 



(<r • %) <r, (1) 



vj + — — (ct ■ Vij) d-, (2) 



where = v^ — Vj is the relative velocity and a is the unit vector joining the center of the 
two particles at contact. The z axis will be taken in the direction of the gravitational field, 
so that the particles are submitted to an external force of the form / = —mgoe z , go being 
a positive constant and e z the positive unit vector along the z axis. Collisions of particles 
with the movable piston on the top are also considered as smooth and inelastic, otp being 
the coefficient of normal restitution for them. Therefore, in a collision between particle i 
and the piston, the component v i>x of the velocity of the particle perpendicular to the z axis 
remains unchanged, 

Vi,x = V ^ ( 3 ) 

while the component Vi )Z of the velocity of the particle and the velocity V z of the piston are 
instantaneously modified accordingly with 

<* = v ^ t 1 + ap ^ v ^ ~ V ^ ^ 

m + M 



and 



TYl 

m + M 



respectively. When the transversal section of the system, i.e. the size of the piston W, 



is smaller than a critical value 
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16], the gas reaches, after a transient time interval, a 
stationary state with gradients only in the direction of the gravitational field. If, in addition, 
the inelasticity of the system is small, the inelastic hydrodynamic Navier-Stokes equations 



with the appropriate boundary condition, M, M P-ide an aeenrate description of the 
stationary state [19) . As the coefficient of normal restitution of the gas a decreases, signif- 
icant deviations from the predictions following from the Navier-Stokes equations show up. 



They are due to the coupling between inelasticity and gradients that exists in the stationary 
state, in such a way that strong inelasticity implies large gradients of the hydrodynamic 
fields. This coupling is a peculiarity of the steady states of inelastic fluids, following from 
the balance between the energy dissipated because of inelastic cooling and the hydrodynamic 
energy fluxes. 

A previous analysis, carried out in jl9], focussed on the macroscopic description of the 
granular gas in terms of the density and granular temperature fields, the velocity field being 
zero. At this level of description, the role of the movable piston on top of the gas is to partially 
determine the boundary conditions needed to solve the hydrodynamic equations for the 
steady state under consideration. Here, the interest will be on the fluctuations of the movable 
piston, namely on its position fluctuations. Some results for the velocity fluctuations have 
been reported elsewhere [jjj]. There, it was shown that the steady state velocity fluctuations 
of the piston are gaussian with zero mean for otp> 0.6 and a > 0.8. Nevertheless, no simple 
relationship between the second moments of the velocity distributions of the piston and the 
gas next to it was found. It is worth to remark that there is no reason to expect such a 
relation to exist at a macroscopic level of description, i.e. involving only the hydrodynamic 
fields and the parameters of the system. Actually, the simulation results reported in [ijj 
indicate that the details of the velocity distribution of the gas, beyond its first few moments, 
are relevant to determine the second moment of the velocity distribution of the piston. 

It is clear that the position fluctuations of the movable piston are related with the volume 
fluctuations of the inelastic gas. Actually, this relationship can be made direct and exact 
by properly choosing the nature of the vibrating wall located at the bottom of the system. 
The mission of the latter is to energize the system, keeping the particles fluidized. The 
expectation is that the behavior in the bulk of the system is independent of the details of 
the way in which this wall is being vibrated. Consequently, the simplest possible choice has 
been used in all the results to be reported in the following. The bottom wall is vibrated with 
a sawtooth velocity profile, having a velocity Vw- This means that all the particles colliding 
with the wall find it moving upwards with that velocity 20|, |2l|. Besides, the amplitude 
of the wall motion is considered much smaller than the mean free path of the particles in 
its vicinity, so that the position of the wall can be taken as fixed at z = 0. Therefore, the 
dynamics of the vibrating wall at the bottom does not induce directly any change in the 
volume (or area) occupied by the granular gas. Also, and again for the sake of simplicity, 
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collisions of the particles with this wall will be considered as elastic. 

In the event-driven simulations carried out, periodic boundary conditions were used in the 
x direction. The width of the system and the number of particles in it were fixed to W = 70d 
and N z = N/W = 6<i _1 , respectively. Moreover, the velocity of the vibrating wall vw was 
chosen in each case large enough, not only to fluidize the system, but also to guarantee 
that the density remains small throughout the granular gas and, consequently, the dilute 
limit can be expected to be accurate. In this case, the dependence of the hydrodynamic 
profiles on vw is very simple and follows by dimensional analysis 19|. The dependence of 
the position fluctuations of the piston on this velocity will be discussed later on. The value 
of the coefficient of normal restitution for the collisions between particles has been varied 
within the interval 0.85 < a < 1. This includes a range of values for which the Navier- 
Stokes hydrodynamic description is not accurate, due to the coupling between inelasticity 
and gradients already pointed out [19]. The coefficient of restitution for the collisions of the 
particles with the movable piston has been set to «p = 0.99 in most of the simulations being 
reported, but it has been verified that the results depend very weakly on the value of this 
coefficient, remaining practically the same when ap is decreased, at least up to ap = 0.8. 
Some examples of this behavior will be given below. 

In all the simulations being presented, it was observed that the height Z of the movable 
piston oscillates about an average value < Z >= L, once the steady state is reached. As 
an example, in Fig. [T]the time evolution of the scaled position of the piston, Z* = Zgo/v^, 
is shown for a = 0.95, ap = 0.99, and three choices of the mass of the piston: M = 30m, 
M = 75m, and M = 150m. Time r is measured in accumulated number of collisions per 
particle. It is observed that both, the value of L and the amplitude of the fluctuations, 
decrease as M increases. Of course, this is the expected behavior. Also notice that the 
trajectory of the piston does not exhibit systematic oscillations, but it is apparently random. 
This indicates that the motion of the piston does not have any hydrodynamic component 
induced, for instance, by the vibrating wall at the bottom. 

From the steady trajectory of the piston, the probability distribution for its position can 
be built. To increase the statistics, several trajectories have been generated for each set of 
values of the parameters. As already mentioned, the fluctuations of the piston depend on 
the parameters defining the system. To see whether this dependence occurs only through 
the first two moments of the probability distribution, a normalized length i has been defined 
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FIG. 1: Time evolution of the dimensionless scaled position Z* = Zgo/vyr of the movable piston 
located on top of the system, once the steady state has been reached. There are 420 disks in a 
box of width W = 70d. In all cases the coefficient of normal restitution of the particle collisions 
is a = 0.95 and that for the particle-piston collisions is ap = 0.99. Results for three values of 
the mass of the piston M are shown, as indicated. Time t is measured in accumulated number of 
collisions per particle. 



as 



Z - L Z* - L* 



(6) 

°z z 

where oz is the square root of the second central moment or standard deviation of Z , a\ =< 
Z 2 > —L 2 , and the star indicates that lengths are being measured in the dimensionless scale 
defined above. In Fig. |2l the obtained probability distribution of £, P(£), is plotted for a 
system with a = 0.94. Again, three values of the mass ratio have been considered, namely 
M = 36m, M = 60m, and M = 120m. It is seen that the probability distributions are 
accurately fitted by a Gaussian (solid lines), at least up to values of the probability density 
of the order of 10 -4 . A similar behavior has been found in all the simulated systems with 
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TABLE I: Third and fourth moments of the position distribution of the piston, obtained from the 



simulations, ine position L 


' is given in the dimensionless scale defined in Eq. ([6]). 
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parameters within the ranges mentioned above, although it seems that a small but systematic 
deviation shows up as the mass of the piston M becomes smaller, approaching the mass of 
the particles. A possible explanation for this behavior is that, as the mass of the piston 
decreases, the amplitude of its position fluctuations increases, and the effect of the external 
gravitational field breaks the symmetry of the fluctuations around the average position. To 
check this idea and to quantify the deviations from the Gaussian of the position fluctuations, 
the third and fourth moments of i have been computed from the simulation data. The 
results for some of the simulations are given in Table [B For a Gaussian distribution it is 
< £ 3 >= and < £ 4 >= 3. The deviations of the third moment from the Gaussian value 
are much stronger than those of the fourth one, supporting the idea that the main cause of 
the deviation from the Gaussian is due to the symmetry breaking produced by the external 
field, when the mass of the piston is not much larger than the mass of the particles. In any 
case, the deviations are rather weak and it can be concluded that, in the explored parameter 
region, the position fluctuations of the piston can be considered as Gaussian with a very 
good accuracy 
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FIG. 2: Steady position distribution of the piston P(l) in both normal and logarithmic scales. The 
symbols are from the simulations while the solid lines are Gaussian with unity dispersion. The 
data correspond to three systems differing in the mass of the piston, as indicated. The coefficient 
of normal restitution for the particle collisions is a = 0.94. The dimensionless length £ is defined 



19| it was shown that the values of the average position of the piston L scale with 



In ref . 

v w in the low density limit and, consequently, L* does not depend on the velocity of the 
vibrating wall in this limit. The extension of this result to volume fluctuations requires to 
go beyond hydrodynamics. To investigate whether <j* z also has the above scaling property, 
several series of simulations have been performed varying the value of v\y, while keeping 
constant all the other parameters. It is important to stress that Vw was always chosen large 
enough as to fluidize the system and to avoid the presence of regions with density above what 
is considered the low density range. In all the cases investigated, there was no dependence 
of o* z on v-\y, within the statistical uncertainties. Therefore, in the low density limit, the 
standard deviation of the piston position seems to scale with the square of the velocity of 
the vibrating wall, i.e. in the same way as the average position L. 

In Fig. [31 the relative standard deviation az/L = o* z /L* is plotted as a function of the 
mass ratio M/m, for several values of the coefficient of normal restitution of the gas a in the 
interval 0.85 < a < 0.98. The coefficient of normal restitution for particle-movable piston 
collisions is in all the cases ap = 0.99. It is observed that, for given a, there is a region 
in which az/L decreases as the mass ratio increases. This effect is less pronounced the 
smaller the coefficient of restitution a, i.e. the more inelastic the collisions. For large values 
of M/m, the results in the figure indicate that az/L tends to a plateau with a constant 
value. The values of the mass ratio needed to reach the plateau monotonically decrease 
as a decreases. The data in Fig. [3] also indicate that for constant mass ratio, the relative 
fluctuations increase as a decreases. 

To show that the influence of the inelasticity of the particle-movable piston collisions is 
much weaker than that of the particle-particle collisions, in Fig. H] the relative standard 
deviation is plotted as a function of the mass ratio for two pairs of data. Each pair corre- 
sponds to the same value of a (namely, 0.98 and 0.94), but different values of ap (namely, 
0.8 and 0.99). Although the variation of ap is almost five times the variation of a, it is 
seen that the data corresponding to the same a are much closer than those corresponding 
to the same ap. On the other hand, there is a relevant qualitative feature to be stressed. 
While decreasing a produces an increase of the relative fluctuations, decreasing ap has the 
opposite effect: the relative fluctuations also decrease. 
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FIG. 3: (Color online) Dimensionless relative standard deviation crz/L of the position of the 
movable piston as a function of the mass ratio M/m, for several values of the coefficient of normal 
restitution of the gas a, as indicated in the insert. The curves are guides for the eye. In all the 
cases, the coefficient of restitution for the gas-movable piston collisions is ap = 0.99. 

III. COMPRESSIBILITY AND EFFECTIVE TEMPERATURE 

To measure the facility of the system to be compressed, define a coefficient of compress- 
ibility k by 

where pi is the pressure of the granular gas in the vicinity of the movable piston and V is the 
volume (area) of the system. The derivative in the above equation is computed at constant 
value of all the parameters denning the system, a, ap, N z and Vw, except pi = Mg /W, 
as it follows from the definition of the pressure. In the following, pi will be modified by 
changing the mass of the piston M, keeping go and W unchanged. The reason for this 
choice is twofold. First, changing g is equivalent to modifying vw and, second, increasing 
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FIG. 4: (Color online) Dimensionless relative standard deviation crz/L of the position of the 
movable piston as a function of the mass ratio M/m. The symbols are simulation data and the 
lines guides for the eye. The two upper curves correspond to a = 0.94, and the two lower ones to 
a = 0.98. In each case, two values of ap have been employed: 0.99 (circles) and 0.8 (triangles). 

W can lead to the set up of transversal instabilities, as already mentioned. The notation 
stresses the constancy of the velocity of the vibrating wall. Although other compressibility 
coefficients could be denned, the one in ([7]) has the advantage of being easy to implement in 
experiments, at least at a conceptual level. On the other hand, it is worth to remark that 
by keeping vw constant and changing M, both the hydrodynamic profiles inside the fluid 
and the power injected into it through the vibrating wall are modified. In particular, the 
latter is given by [3] 

Qo = W (n z + -^Lj mg v w . (8) 

Therefore, increasing M while keeping vy/ constant produces an increase of the power Q . 
For the two dimensional systems being considered here, the definition (J7J) is equivalent to 
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The idea of introducing a compressibility for vibrated granular fluid was already used in 



ref. 



22] in the context of hydrodynamical stability analysis. In a molecular system at 



equilibrium, the isothermal compressibility 



< V > V dp 
is related with the volume fluctuations by 



* = -^f^l do) 

T 



a\ =< {V- < V >) 2 >= k B T < V > k T , (11) 

where kp is the Boltzmann constant. In non-equilibrium states, there is no reason to expect 
the above relationship to hold, but it is tempting to employ it to define an effective temper- 
ature parameter of the system, T e ff, expecting it to have some intrinsic physical meaning. 
Therefore, taking into account that the volume fluctuations in the case being considered are 
associated to fluctuations of the position Z of the movable piston, we define T e ff through 

2 T eff ( dL 



Tlsff^- (12) 

The Boltzmann constant has been set equal to unity as it is usually done when defining 
the granular temperature from the average kinetic energy of the grains. From Eq. (1I2"|) it 
follows that T e ff relates the volume response to a pressure perturbations with the steady 
volume fluctuations of the system. Of course, Eq. (TT2]) by itself is just a mathematical 
definition and does not add anything to the physical understanding of the system. On the 
other hand, the definition would become relevant if this effective temperature were related 
to other temperature-like parameters of the system. Two main candidates clearly stand out: 
the granular temperature of the gas in the vicinity of the piston , TV , and the temperature 
parameter of the piston, Tp, defined through Tp = M < V 2 > [19]. Both parameters, Tl 
and Tp, are not at all the same, as it should be the case if the system under consideration 
were at equilibrium and energy equipartition would apply. Violation of equipartition is a 



general feature of granular systems known since long ago 



23J and that has attracted a lot 



of attention in the last years. A more detailed discussion of this issue for the set up being 
considered here is given in jl4| . In Fig. [5j event driven simulation results for the ratio Tp/T^ 
are plotted as a function of M/m. Several values of the coefficient of normal restitution of 
the gas in the interval 0.85 < a < 0.98 have been considered, while again ap = 0.99 for all 
the data shown. It is observed that the behavior of this temperature ratio is quite intricate. 
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FIG. 5: (Color online) Ratio of the temperature parameter of the piston, Tp, to the granular 
temperature of the gas next to it, Tl, versus the mass ratio M/m, for several values of the coefficient 
of restitution of the gas a, as indicated in the insert. The curves are guides for the eye. In all the 
cases, the coefficient of restitution for the gas-movable piston collisions is ap = 0.99. 

For instance, Tp/T^ decreases as M/m increases for a > 0.95, but it happens the other 
way round for a < 0.95. Actually, as a decreases below this value, the increase of the 
temperature ratio is rather fast, and Tp reaches to be up to four times larger than Tl- If the 
parameters Tp and Tl were interpreted as real temperature parameters, this latter behavior 
would be fully counterintuitive. The temperature of the heated body (the piston) is larger 
than the temperature of the heating one (the gas next to the piston). 

From the values of oz and (dL/dM) Vw obtained from the event driven simulation data, 
the effective temperature T e ff has been computed by means of its definition, Eq. (|12p . In 
the Appendix some details are given of the way in which the above derivative was actually 
evaluated. Then, in Figs. |6] and [7] the temperature ratios T e ff/Tt and T e ff/Tp are plotted, 
respectively, as a function of the mass ratio for the same systems as in Fig. [5j A clear 
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FIG. 6: (Color online) Ratio of the effective temperature, T e jf, denned in Eq. (|12|) to the tem- 
perature of the gas in the vicinity of the piston, Tl, as a function of the mass ratio for the same 
systems as in Fig. 

difference is observed in the behavior of the two temperature ratios for a given constant 
coefficient of restitution a. While T e ff/TL exhibits a strong dependence on M/m and does 
not seem to tend to a well denned limit as it increases, the dependence of T e ff/Tp on 
M/m is very weak, being only appreciable for the least inelastic cases and when the mass 
ratio is small. Thus it is concluded that for large enough mass ratio M/m, the effective 
temperature is proportional to the temperature parameter of the piston with a coefficient 
of proportionality that is independent of the mass of the piston, i.e. 

T eff = b(a,a P )T P , (13) 

where a possible dependence of the coefficient on ap has been included. 

To identify the dependence of the coefficient b on a, in Fig. [H] the plateau values of 
T e ff/Tp, reached upon increasing the value of the mass ratio M/m, are plotted versus 
1 — a 2 . The points are very well fitted by a straight line, indicating that T e fj/T P grows 
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FIG. 7: (Color online) Ratio of the effective temperature, T e fj, defined in Eq. (fT2j) to the temper- 
ature parameter of the piston, Tp, as a function of the mass ratio for the same systems as in Fig. 
5. 

linearly with 1 — a 2 , at least in the considered interval, 0.85 < a < 0.98. It is worth to note 
that some care is needed when extrapolating to the elastic limit a — > 1 these results. In 
this limit, a stationary state is only possible if, in addition, the vibrating wall is arrested, 
i.e. also the limit vw — )• is taken. But all the previous discussion has been carried out at 
constant velocity of the vibrating wall. This explains why extrapolation of the linear fitting 
in Fig. E] does not lead to b = 1, and this does not mean any kind of contradiction. 

In Fig. the dimensionless compressibility, —mL~ l (dL/dM) vw and the scaled second 
moment of the position fluctuations aerugo/ 'Lb(a,ap)Tp have been plotted as functions of 
M/m. A logarithmic representation is employed. Data for several values of the coefficient a 
have been included. For all the data ap = 0.99. If T e ff had been used instead of b(a, ap)T P , 
the two plotted quantities would be the same by definition, i.e. the filled and empty symbols 
would agree in all the cases. It is seen that the dependence on M/m of the dimensionless 
compressibility is accurately described by a power law of the form (M/m) -3 / 4 , indicated in 
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FIG. 8: Ratio T e ff/Tp for large values of M/m as a function of 1 — a 2 for the same systems as in 
Fig. [5j The symbols are simulation results and the dashed line a linear fit of them. 

the figure by the solid straight line. Nevertheless, it is important to realize that the interval 
of values of M/m for which the above dependence is identified is rather narrow, just one 
order of magnitude, so that its range of validity can be limited. 

IV. DISCUSSION AND SUMMARY 

The aim here has been to investigate the volume fluctuations of a vibrated low density 
gas of inelastic hard disks in presence of gravity, and confined by a movable piston on the 
top. The study has been restricted to the parameter region in which the system reaches 
a steady state with gradients only in the direction of the external gravitational field, i.e. 
perpendicular to the movable piston. In practice, this has limited the values of the co- 
efficient of normal restitution of the gas particles to the interval 0.85 < a < 1. Due to 
the coupling between inelasticity and hydrodynamic gradients, which is peculiar of steady 
states of granular systems, the above limitation also implies restriction to small gradients. 
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FIG. 9: (Color online) Dimensionless compressibility —mL (dL/dM) Vw (filled symbols) and 



scaled position fluctuations of the piston aerugo/ Lb(a, ap)Tp ( empty symbols) as a function of 
the mass ratio, M/m. The different symbols correspond to different values of the coefficient of 
restitution for the collision between particles, a, as indicated in the insert. The straight line has a 
slope —3/4, and it is a guide for the eye. 



Nevertheless, the analysis carried out in ref. 



191 ] indicates that the range of hydrodynamic 



gradients considered here exceeds the limit of validity of the Navier-Stokes approximation. 

Some of the main results can be summarized as follows: i) for large mass of the movable 
piston compared with the mass of the gas particles, the volume fluctuations are Gaussian 
with very good accuracy, ii) the square root of the second moment of their distribution 
scales with the square of the velocity of the vibrating wall at the bottom, i.e. in the same 
way as the amplitudes of the hydrodynamic fields in the gas, iii) by requiring the same 
relation between volume fluctuations and compressibility as in equilibrium systems to be 
verified, an "effective temperature" can be defined, iv) the effective temperature turns out 
to be proportional to the second moment of the velocity fluctuations of the piston, with 
a proportionality parameter that depends on the inelasticity of both the particle-particle 
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and particle-piston collisions, but it seems to be independent of the mass of the piston, and 
v) the effective temperature can not be related in a simple way to the temperature of the 
granular gas; even more, the relationship between both parameters is not monotonic. 

A relevant open question is the relationship between the granular temperature of the gas 
in the vicinity of the piston and the temperature parameter of the piston, the latter defined 
from the second moment of its velocity distribution. An explanation of the simulation results 
seems to require a detailed knowledge of the velocity distribution function of the gas next 
to the piston [fjj]. If this is the case, approximated solutions of the Boltzmann equation, as 
provided by instance by the Chapman-Enskog procedure in the first Sonine approximation, 
would not be of enough accuracy as to describe the deviation from equipartition between 
the gas and the movable piston. 

The present study complements the one in ref. [l^], in which the velocity fluctuations of 
the piston were investigated in detail. A natural issue now is whether the velocity fluctu- 
ations and the position fluctuations of the piston are correlated. We have computed from 
the simulation data the joint probability distribution for the position and velocity of the 
piston and compared it with the product of the marginal distributions for the position and 
the velocity. Both results agree within the statistical uncertainties, indicating the absence 
of correlations. 
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Appendix A 

In ref. 



19], an expression for the average position of the piston, L, was derived by using 
a hydrodynamic description of the granular with the appropriate boundary conditions. The 
theoretical prediction was showed to be in reasonable agreement with the simulation results. 
Although this expression could have been used to compute (dL/dM) Vw , here this quantity 
has been obtained from the simulation data for the sake of consistency. 

When trying to compute from the values of L as a function of M/m the derivative of 
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the former with respect the latter, the technical problem arises that the considered values 
of M/m are separated by intervals of the order of 10 or 15 units. Decreasing this interval 
would require to strongly increase the simulation time. In addition, some smoothing process 
should be used. Here, a different approach has been followed. The simulation data for L have 
been fitted by an analytical function of the mass ratio whose derivative has been afterwards 
computed. The trivial choice of a polynomial in m/M does not work so, motivated by the 
equilibrium elastic result, an expression of the form 

mv 2 w V M J 

where A and B are two adjustable dimensionless parameters, was used. It turned out that 
it fits very well the results for all the values of a and ap considered here. Moreover, the 
values of the two fitting parameters are quite stable, in the sense that their value seem to 
converge very fast when the number of fitted data is increased. 
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